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Cn| Abstract 

^ We consider the refinement of the holographic entanglement entropy on a disk region for the 



holographic dual theories to the AdS solitons and AdS black holes, including the corrected ones 
by the Gauss-Bonnet term. The refinement is obtained by extracting the UV-independent piece of 
the holographic entanglement entropy, the so-called renormalized entanglement entropy which is 



o independent of the choices of UV cutoff. Our main results are (i) the renormalized entanglement 

,^ entropy of the AdS^^+i soliton decreases monotonically for d = 3, however for d = 4, 5 violation 

_^ of the C-theorem occurs around the deconfinement /confinement phase transition; (ii) there is no 

> 

^y-^ topological entanglement entropy for AdSs soliton even with Gauss-Bonnet correction; (iii) for 

in 

T— I the AdS black holes, the renormalized entanglement entropy obeys an expected volume law at IR 

-yJ regime, and the transition between UV and IR regimes is a smooth crossover even with Gauss- 

pq Bonnet correction; (iv) based on AdS/MERA conjecture, we postulate that the IR fixed-point state 



for the non-extremal AdS soliton is a trivial product state. 
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I. INTRODUCTION 

Quantum entanglement is an important theoretical probe to understand some particular 
feature of the strongly coupled systems [U |2] , such as the topological ordered phases which 
are believed to be related to the long-range entanglement [31 1^. On the other hand, the 
nature of short-range entanglement for generic ground states yields the famous area law 
[5l E] . The entanglement entropy is plagued by the UV cutoff, however, some of the encoded 



information is related to the counting of number of degrees of freedom. The famous example 
is the entanglement entropy of the (1 + l)-dimensional conformal field theory (CFT), for 
which the coefficient of the logarithmic UV divergent term is proportional to the central 
charge of the CFT. Thus, one task for the physical interpretation of the entanglement 
entropy is to extract such kind of the UV-independent piece, or so-called the renormalized 
entanglement entropy. If these renormalized entanglement entropies are indeed related to 
the number of the effective degrees of freedom [T71 |20l |2T], then one may expect that they 
should obey some sort of C- or F-theorem along the renormalization group (RG) flow, at 
least for relativistic quantum fleld theories [2SH2B]- 

However, it is difficult to evaluate the entanglement entropy directly even in the text of 
free fleld theory, which is usually based on replica method [H [8] , not mentioning to evaluate 
it directly for the strongly coupled theory. Fortunately, it was proposed in [9HTT] that 
in the context of AdS/CFT correspondence, the holographic entanglement entropy has a 
simple geometric representation, which is the area of the minimal hyper surface in the bulk 
with its UV boundary coincident with the entangling surface in the dual fleld theory. As 
usual, the holographic entanglement entropy is plagued by the UV cutoff, and one should be 
careful to extract the UV-independent piece which is free of the UV cutoff ambiguity. The 
explicit calculation of the entanglement entropy of the d-dimensional relativistic conformal 
fleld theory (CFT) with its holographic dual the gravity in ((i + l)-dimensional anti-de Sitter 
(AdSd+i) space, gives the following generic UV scaling structures [TT] 

te + --- + f + const + ;|---, dodd, 

[^ + ■■■ + $ + const log f + ;J ■ ■ ■ , d even, 

where R is the linear size of the entangling surface (in this paper we will consider the disk 
case), and e is the UV cutoff. This UV structure is consistent with the one obtained from 
the fact that the entanglement entropy should be an even function of extrinsic curvature of 
the entangling surface |19]. Moreover, the constant parts in the above are UV-independent 
pieces, which will not change under the redeflnition of the UV cutoff e, and should be 
identifled as the renormalized entanglement entropies. 

Recently, it is proposed in [181 ES] how to extract from ([I]) the renormalized entanglement 
entropy. The basic idea is to construct some rf-dependent function fd^Rdn) of differential 
operator Rdji so that when acting on ([I]) by this operator one will extract the aforementioned 



UV-independent pieces, which are related to the central charges of the CFTs and should 
be positive. The detailed form of fa is given in [22] • One then applies the same differential 
operator /^ to the entanglement entropy of the non-CFTs and extracts the corresponding UV 
independent pieces, which should be the C-functions and are expected to be monotonically 
decreasing along the RG flow as R increases. Similar works have recently been done in 
[23] 121] . However, the way of extracting the renormalized entanglement entropy is far from 
unique as in the usual case for other renormalized quantities plagued by UV divergence. 
Despite that, for the extracting quantities to be related to the number of the underlying 
degrees of freedom, we should require it to be positive and obey some C-theorem at least 
at the very beginning of RG flow. The aforementioned fd^RdRj is devised to satisfy these 
constraints, and is succinct and scale-adaptive. 

In this paper, we would like to generalize the above extraction scheme to the one for the 
holographic dual non-CFTs which are gapped (and non-relativisitc) or finite temperature 
version of CFTs (also non-relativisitic due to the temperature and chemical potential), and 
check if the resultant renormalized entanglement entropies obey the C-theorem or not. Our 
motivation for considering such cases is partly to see if the extraction scheme is universal or 
not. On the other hand, there may have topological order for gapped systems, which can 
be encoded in the constant piece of the renormalized entanglement entropy, the so-called 
topological entanglement entropy |3l HI [19]. We will like to examine its existence by the 
aforementioned extraction scheme for the holographic duals considered here q 

The holographic duals of the finite temperature version of CFTs are the black holes in 
AdS spaces, the UV divergence structure of which is similar to ([I]), thus we can apply the 
same differential operator given in [22] denoted by /^ iRdn) to extract the renormalized 
entanglement entropy. As we shall see, the renormalized entanglement entropy shows an 
expected smooth crossover from the UV regime to the volume law in the IR regime [35j , the 
latter captures the extensiveness of the thermal entropy encoded by the black hole horizon 
as the entangling surface. 

On the other hand, the holographic duals of the gapped systems considered here are the 
so-called AdS solitons, which can be obtained by double Wick rotation of the AdS black 
hole metric and then by compactifying one of the transverse dimensions. The warped size 
^ For earlier studies on the UV structure of the holographic entanglement entropy for the AdS solitons, see 

[IMS] 



of the compact circle shrinks to zero at some finite value of AdS radial coordinate so that it 
caps out the rest of the original AdS geometry. The capped geometry implies an IR fixed 
point of the dual deformed CFT at finite energy scale by the UV-IR correspondence, thus 
it is dual to a gapped system. Moreover, this compact dimension is a spectator for the dual 
deformed CFT, i.e., the entangling surface wraps over it, it then yields different UV scaling 
structure from ([I]). Instead, for AdS^+i soliton it looks like 

where Lg is the fixed proper size of the compact circle, and S\jy is the UV structure of AdS^+i 
given in ([I]). We shall then adopt a differential operator to extract the UV-independent piece 
of the entanglement entropy, denoted as S'uv-ind, also called the "renormalized entanglement 
entropy" for short. Since the two UV scaling structures are related, it is straightforward to 
see that the differential operator 

g,{Rdn) := ^ff''\Rdn)R (3) 

will retain the salient feature of /^ proposed in [22], namely, being succinct and scale- 



U 
adaptive, and to result in a positive C-function obeying C-theorem at the UV regime. For 

example, in the UV limit, the renormalized entanglement entropy extracted from AdS^^+i 

soliton is q/-R where q is the renormalized entanglement entropy extracted from pure 

AdSrf+i space (central charge of the dual CFT). After some manipulations, we can write 

down the RG flow of the renormalized entanglement entropy as follows 

"^^ 1 (^ (RjR + mii ■■■{R^-{d- 4)) i, d even, 

where 5* is the holographic entanglement entropy for the AdS^+i soliton, whose UV scaling 
behavior is related to the one for AdS^+i space via (J2]). Despite the nice feature in the 
UV-regime, we may not expect the C-theorem to be held along the RG flow due to the 
non-relativistic nature of the gapped systems. This makes the evaluation of the RG flow 
of the renormalized entanglement entropy for the gapped systems an interesting task. It is 
indeed the main goal of this paper. 

Besides, if there is no topological order, one may expect the IR flxed point of the gapped 
system will be a trivial product state after performing proper local unitary transformation 




FIG. 1. Upper: The procedure of MERA or equivalently quantum state renormalization group 
transformation for the gapped system. The circle at each step denotes the surface enclosing the 
chosen region, and the links crossing it denote the entangled pairs which contribute to the entan- 
glement entropy after tracing out the wave function outside/inside the circle. The length of the 
link is the distance between the entangled pair, and signifies the entanglement at that length scale. 
At each step of MERA, the entanglements at the corresponding scale are removed. There are two 
possible end states at the IR fixed-point: (4a) the trivial product state and (4b) the entangled 
state protected by symmetry or topological order. Lower: The corresponding holographic minimal 
surfaces in the bulk AdS soliton. The (4a) and (4b) in MERA yield the minimal surfaces of cylin- 
der and disk topologies, respectively. Moreover, the entanglement entropy at each scale of MERA 
is encoded in the area of the minimal surface above the yellow bar at that scale. As seen, such 
area for (4a) is negligible compared to (4b). It then suggests that (4a) is a product state without 
entanglement but (4b) is not. More detailed explanation will be given in section 5. 



to remove the short range entanglement. If so, it imphes that the rate of change of the renor- 
mahzaed entanglement entropy along the RG flow approaches to zero in the IR limit. One 
may also reveal this kind of feature geometrically in the holographic dual gravity. We will 
try to argue this is indeed the case based on the proposal of AdS/MERA (multi-scale entan- 
glement renormalization ansatz) [50l [51] by just looking into the dominant topology of the 
large holographic entangling hypersurfaces. We briefly summarize the idea of AdS/MERA 
and the associated entangled nature of IR fixed-point state in Fig. [T| and the more detailed 
explanation will be given in section |Vj 

Our paper is organized as follows. In section [TT], we will extract the UV-independent piece 
of the holographic entanglement entropy for the AdS^+i soliton with generic form of metrics. 
Then, we will evaluate numerically the RG behavior of the UV-independent piece. We also 
discuss how to extract the topological entanglement entropy from the UV-independent piece. 



In section III, a similar consideration goes for AdS black holes. In section [TV] we will extract 
the UV-independent piece of the entanglement entropy and its RG flow for the AdSs soliton 
and black hole corrected by the Gauss-Bonnet term. We then conclude our paper in section 
\V\ by discussing the entangling nature of the IR fixed-point state of the holographic dual 
theory based on the proposal of AdS/MERA. 

II. HOLOGRAPHIC RENORMALIZED ENTANGLEMENT ENTROPY FOR ADS 
SOLITONS 

In this section, we will first discuss how to extract the UV-independent piece of the 
entanglement entropy for the AdS soliton, which is free of the UV cutoff and the associated 
ambiguity. Then we will discuss how to extract the topological entanglement entropy from 
the UV-independent piece, which should be encoded in the constant piece in its IR limit. 

We will consider the AdS soliton with following form of metrics in the Poincare coordi- 
nates, which can be obtained from the double Wick rotation of some asymptotically AdS 
space: 

ds' = ^ [^ + f{z)de' - de + dr^ + rHn,^ , (5) 

where the harmonic function f{z) can take the general form as follows 



'/ n=l 



f{z) = [i-kA\\i-h^^]{i + y] c^z-). (6) 



We assume the c„'s are chosen appropriately such that 1 + J2n=i ^nZ^ does not contain poles 
and zeros aX z = zq. The parameters ki and ^2 can be tuned to yield different IR behaviors. 
The metrics include the pure AdS space by choosing A;i = ^2 = and setting zq = 1. 

The simplest AdS soliton is the one with ki = 1 and k2 = —1 and with c„ chosen so that 
f{z) = 1 — (— )^^'^. By choosing the proper period of ^-coordinate, denoted by Lg to remove 
the conical singularity, this metric has a smooth tip at z = Zq which corresponds to the IR 
gap of the dual theory. Note that ggg is different from the other components rather than g^z- 
This implies that the boundary space-time, which includes also ^-coordinate, is not Lorentz 
invariant. Moreover, the proper size ^/gggLg of the ^-direction depends on the RG scale z 
so that it yields a d-dimensional UV theory but a (rf — l)-dimensional IR theory since the 
proper size of 6 shrinks to zero there. One can also turn on some deformation operators 
to the dual boundary theory of pure AdS soliton, which are encoded in c„'s capturing the 
deviation from the ones for f{z) = 1 — (— )^~'^. For example, one can double Wick rotate 
the AdSs charged black hole with the harmonic function f{z) = 1 — mz^ + q'^z^. This is then 
dual to a boundary theory with non-zero current density condensate or magnetic fluxes. 
More complicated case can be obtained from other deformations of the pure AdS metric, 
such as the hairy scalar AdS black hole [36| |37] or even AdS R-charged black hole [30] . 

For simplicity, we will set LacIS = 1 ^^^ focus on d = 4 and d = 5 case, but also including 
d = 3 case for completeness. Here we refer d to the space-time dimension of the UV theory. 
In some literature, it refers instead to the space-time dimension of the IR theory, which is 
one dimension less than the UV one. 



A. Extracting the renormalized entanglement entropy 

To evaluate the holographic entanglement entropy, one should find out the minimal sur- 
face with its boundary enclosing the entangling surface. This is done by finding the solution 
of the equation of motion derived from the action for the area of the above hyper- surface, 
i.e., 

v/det(7ind = ^d-3 / dz-^,/l + fr^:=nd-3 dz C, (7) 

8 



where (^ind is the induced metric on the hyper-surface, and t^ = %■ The holographic entan- 



glement entropy S is related to the area Ahy S = -^^ 
The equation of motion for r(z) explicitly is 



A 



■0 



1{d - 1)/Vr'' + 1z{d - 3 - rff) + /(2(ci - 'i)zf' - r{-2{d - l)r + zff' + 2zr)) = 0, (8) 

where / = ^^. The minimal surface will have different IR behaviors depending on the 
linear size R. For generic AdS soliton metric, the small R minimal surface will have a disk 
topology and z^ is the turning point such that r{zm) = 0. On the other hand, the large R 
one will end on the z = zq, thus Zm = zq with a cylinder topology, see Fig.[2j However, for 
the case with extremal harmonic function, i.e., ki = k2 = 1, only exists disk topology for all 
R. 
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FIG. 2. Disk (blue) and cylinder (red) topology of the minimal surface for AdS soliton. 



Varying A with respect to R with z = e fixed, and using the Hamilton- Jacobi method, 
we find that |22] 






Hie 



drje) 
dR 



-H(e 



dr{e) 
dR 



where 



H:= 



5C^ 

6r 



/vlOi CJ -U ly. 



^d-3 



Z 



d-1 



v/r+T^ 



n = uf-c 



'^-lyiTT^' 



z 



(9) 



(10) 



^ For simplicity, hereafter we will omit the angular factor fld-3 and will not distinguish between A and 
A/fld-3 and similarly for the quantities related to A such as S, S'finitc and S'uv-ind- 



The first term in the first equahty of ^ is dropped because of the IR boundary condition 
for the minimal surface, i.e., 

r{ztn) = s.t. l-L^Zm) = for disk topology, (11) 

dzm dzo 



dR dR 



for cylinder topology. (12) 



Note that ^ only depends on the UV behavior of the solution r(z). So the resulting scaling 
behavior should hold for both disk and cylinder topologies. However, since the UV boundary 
condition alone cannot determine the full solution, some IR information will be encoded in 
r{z) and affect the renormalized entanglement entropy implicitly. 

Therefore, we only need to extract the UV behavior of the solution r{z) to yield ^, from 
which we can obtain the RG flow of the holographic entanglement entropy after subtracting 
off the UV divergence and its associated ambiguity. We postulate the UV behavior of the 
solution r{z) as 

r{z) = i? + 6o log I + J](a„ + 6„ log ^)z^ , (13) 

n=l 

we then plug it into (Isj) to determine a^'s and 6„'s. 



1. AdS^ soliton 



For concreteness, we consider d = 4 case first. We find that 

where ■ ■ ■ denotes the higher order terms which can be determined by 04, /cj's and c„'s but 
are not relevant for our purpose. An important point is that the equation of motion at the 
UV expansion can not determine 04 (i?). Instead one should determine it by solving the 
full equation of motion. In other word, a^lR) encodes some IR information of the minimal 
surface and the nontrivial RG flow of the holographic entanglement entropy. Especially, it 
should tell when the phase transition occurs between disk and cylinder topologies by tuning 
R. This phase transition is nothing but the deconfinment /confinement phase transition [12] 
with disk topology corresponding to deconfined phase at small R, and the cylinder one to 
the confined phase at large R [131 E] ■ 

10 



Plugging (14) into ^, we obtain 
dA _ ^ ^ ^ ^^ ^ -kf -kl- kik2 + Ci{ki + k2) -cl + C2 3 

- " ' -2" 

^0 



^^ - -4i?a4(i?) + ^ '- ^^2 ^^ ^-^^+ UV-divergent terms +0{e) 



(15) 
where 0(e) terms vanish at e — t- hmit and are not relevant. The UV-divergent terms are 

which are only defined up to the redefinition of the UV cutoff e. For example, redefining 



e by aoe(l + aie H ) will then shift (16) by terms of C(i?°) and Oi.R-'^) with finite UV 



cutoff-dependent coefficients. This means that the terms of 0{R^) and 0{R~^) in (15) are 
not universal but depend on the UV cutoff. To obtain a UV cutoff-independent refinement 
of the holographic entanglement entropy, i,e, the renormalized entanglement entropy, we 



shall then subtract these kind of terms from (15). 

In this paper, we will consider the differential subtraction scheme given in (pi) and (111). 
To demonstrate how (111) is arrived, we take the current example, i.e., d = A. Using f^^RdR) 
given in [22] and ([3]), after some manipulations the differential operator acting on A to 



extract 5'{fv-ind is 



.(4) 



g^RdR) = \{RdR + l){RdR - 1). (17) 



From 5'uy_;j^jj := g4{Rdji)A and using the commutator relation [Or, ROr] = Or, we can then 
obtain 

'-%^.liW.,2)im./A, (18) 

This is the d = 4 case in (111). 

We numerically solve 04 (i?) for the AdS soliton with f{z) = \ — z"^ and the result is 
shown in the left plot of Fig.[3| in which the blue and red curves denote contributions from 
disk and cylinder topologies, respectively. The 04 (i?) is not single- valued near the phase 
transition between disk and cylinder topology. Since we have no other criterion for picking 
out a preferred value of 04(1?), to remove the additional branches we have to compare the 
on-shell actions of the solutions with both disk and cylinder topologies around the critical 
point. Solutions with the smallest on-shell actions are chosen to be the dominant phase. 

To determine the dominant topology, we introduce S'gnite denoting the finite part of 
the on-shell action. For AdSs solitons, this is obtained by subtracting the divergent part 
•^div ~ 2^ + 8^ i°S ^ numerically from the total on shell action. The S'snite is different 

11 
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FIG. 3. Left: The a4{R) for AdSs soliton with f{z) = 1 — z^. Right: Finite part of on-shell action 
'S'finite for the solutions around the critical point. 



from the S'uv-ind defined in ( 18 ) in the sense that the UV cutoff-dependent terms not being 
removed. In fact, ^finite is related to S'uv-ind via S'uv-ind = fl'4(-R<9ij)Sfinite- On the other 
hand, the S'finite could be used to determine the dominant phase, since it contains the total 
information of the entanglement entropy, up to a divergent part which is the same at every 
value of R for different branches. 

The numerical results of the S'finite are shown in the right plot of Fig. [3lP1 , which indicates 
that for R < 0.703 the disk topology dominates, while for R > 0.703 the cylinder topology 
dominates. Hence the additional branches of ai{R) in the left plot of Fig. [3] are removed. 
The discontinuous jump indicates a quantum phase transition. In fact, we can use the renor- 
malized entanglement entropy Suv-ind as an order parameter to characterize this quantum 
phase transition. 



By using (18) and the results in Fig. 3 we numerically calculate the RG flow of the 



renormalized entanglement entropy. 



dS- 



^ig.Q 



(4) 

dfl '"'' , which is shown in the left plot of Fig.|4l We 

find there is a sharp jump around the critical point, indicating the quantum phase transition 

from the deconfining phase in the UV regime to the confining phase in the IR regime. The 

C-theorem holds in the UV regime, which is expected since we define our subtraction scheme 

in the UV limit. However, the C-theorem is seriously violated around the quantum phase 

transition, indicating an increase of number of degrees of freedom. This is at odds with the 

deconfinement/confinement phase transition. On the other hand, the violation of C-theorem 

^ Note that the S'finite in Fig. [s] is negative such that it cannot be directly interpreted as the entanglement 

entropy or the number of degrees of freedom. Instead, 5uv-ind is positive at least in the UV regime as 

guaranteed by the large initial UV value ^ for small i?, thus it could be thought as the number of degrees 

of freedom. i o 



dSuv-ind/dR 
15 r 
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FIG. 4. Left: The 



dS, 



(4) 



(IR 



for AdSs soliton with f{z) = 1 — z^. Right: The corresponding S^ 



(4) 
UV-ind- 



goes way very rapidly once away from the quantum critical point. In fact, the renormalized 
entanglement entropy remains almost constant in the IR regime, this is consistent with 
expectation for the confining phase or the IR mean field state of a gapped systems, i.e., 
almost all the degrees of freedom are gapped out and the ground state is a trivial product 
state. 

As a check of consistency, we also calculated 5'|jy_;j^j directly through g^iRdRfSfinite and 
the result is shown in the right plot of Fig.|4} It is hard to tell whether the quantum phase 
transition is of first-order or second-order because of the numerical error. The fact that 
'S'uY_j„j is neither monotonic nor positive-definite is similar to the result of GPPZ flow 
obtained in |22] • This might be a generic feature oi d = 4 gapped systems, indicating that 
>S'uv_j„(j might not be a good measure of the number of degree of freedom in four dimensions. 

We now consider the special cases of AdS solitons, the extremal AdS solitons. For such 
cases, there are only solutions with disk topology. Technically, this fact could be realized 
from the IR expansion. Suppose that there exist solutions with cylinder topology, which 
end on z = Zq at Tq = r{zo). We could expand the solution r{z) around r = Tq as following: 



r{z) = ro + di{zo - z) + ^2(^0 - zf + ds^Zo - z)^ + 



(19) 



For non-extremal AdS soliton, we could work out the coefficients di, ^2, d^, ■ ■ ■ order by 
order from the expansion of the equation of motion; however for extremal AdS soliton, one 
find that the coefficients di, d2, d^, ■ ■ ■ turn out to be infinity, which indicates that z'{rQ) 
tends to zero. This means that one can never reach the boundary from z = zq, that is, 
solutions with cylinder topology do not exist. 

On the other hand, the difference from the non-extremal case is that the proper size of 
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the compact circle becomes infinite for tlie extremal one, this means that all the associated 
Kaluza-Klein (KK) modes become massless. That is, the dual field theory is a gapless system, 
and the IR behavior of the disk topology solution reflects this fact. More discussions on this 
will be given in section [V] 
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FIG. 5. Left: The 



dS, 



(4) 



The corresponding 5, 



dR 

(4) 
UV-ind- 



for extremal charged AdSs sohton with f(z) = 1 — 3z^ + 2z^. Right: 



By adopting the same differential subtraction scheme as for the non-extremal case, our 



numerical results for 



dS- 



"^'"'^ as well as 5'{jy_jjjj for extremal charged AdSs soliton with 



f{z) = 1 — 3z^ + 2z^ are shown in Fig. 



:(4) 



Using again the 5'{jy_jjj^ as the order parameter, we 
conclude that there is no phase transition. However, the C-theorem and positive-definiteness 
are still violated for R greater than some specific value, and the IR behavior is also similar 
to the one for the non-extremal case. 



2. AdSe soliton 

Similarly, we now consider the d = 5 AdS soliton. The UV expansion of the solution r{z) 
takes the following form 

z^ 2(ci -ki- k2)z^ 



r(.)=i?-^ + 



+ a4(i?);z* + 05(^)2 + C(2 



where 



04 (-R) 



9Rzo 

-kl - kl - kik2 + ci{ki + k2) - cl + C2 5 



QRzl 



54i?3 



(20) 



(21) 



but a^{R) cannot be determined from the UV expansion and should be solved from the full 
equation of motion. From the above expansion, we obtain 

^ = -bR^a,{R) + 2(ci-fci-/c2) ^ ^2R ^ uV-divergent terms + 0(e), (22) 

CLR oRZq OZg 
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where the coefficient B depends only on the detailed form of the metric, i.e., 

B = cl — 2ciC2+C3 — clki + C2ki + Cikl — kl — clk2+C2k2 + Cikik2 — klk2 + Cikl — kikl — k2- (23) 
The UV-divergent terms take the form 



2R 
3?' 



(24) 
Note that there is no log divergent term as expected for d = odd case. It seems a bit 



miraculous that there is also no 0{l/e) term in (22), however there is such a term if we 



integrate (22) over R. To see this, we substitute (20) into the action ([T]), expand the 



integrand in series of z and then integrate, we will find an additional divergent term 



9e" 



Since it is independent of R, we could not find it in (22). The R scaling behaviors of the 
UV-dependent terms are also different from the relativistic case. 

We numerically solve a^lR) for the AdSg soliton with f{z) = 1 — z^ and the result is 
shown in Fig.[6j Again the blue and red curves denote contributions from disk and cylinder 
topologies, respectively. It is interesting that near the critical point Re ~ 0.9415, there seems 
to be a fractal vortex structure, as is shown on different scales in Fig.[6j This indicates that 
05 (i?) is multi- valued near the critical point. 

To determine the dominant phase, we numerically calculated the finite part of the on-shell 
action with the divergence ^div ~ |^ ~ ^ subtracted, and the result is shown in the right 
plot of Fig.[7| The situation is similar to the right plot of Fig.|3] of the AdSs soliton case, 
though it is a bit hard to distinguish the red and blue curves since they nearly coincide with 
each other. From this plot we read the phase transition point R = 0.9355. For R < 0.9355 
the disk topology is dominant, while for R > 0.9355 the cylinder topology is dominant. The 
05(7?) on the corresponding scale with additional branches removed is shown in the left plot 
of Fig. [7} Note that the fractal vortex structure around Re ~ 0.9415 is totally removed, 
hence it will not bring additional phase transitions. 



d5(5)_. 

The RG flow — ^Jj^""^ is calculated straightforwardly by using (|4|) for d = 5 and is shown 
in the left plot of Fig.|8} The qualitative behavior is similar to the one for AdSs soliton. 
It also indicates the occurrence of the deconfinement /confinement phase transition and the 
violation of the C-theorem in the confining phase. The 5'{jy_;^j is also calculated through 



g5{Rdii)S{inite and is shown in the right plot of Fig.pl which is again neither monotonic nor 
positive-definite, indicating a ffist-order phase transition. 
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FIG. 6. The a5(i?) on different scales for AdSe soliton with f{z) = 1 — z^ : fractal vortex structure 
around the critical point. 
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FIG. 7. Left: The a^{K) in detail for AdSe soliton with /(z) = 1 — z^, with additional branches 
removed. Right: Finite part of on-shell action for solutions around the critical point. 
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FIG. 8. Left: The 



dS, 



(5) 

ST^ for AdSe soliton with f{z) 



1 — z^. Right: The corresponding S, 



(5) 
UV-ind- 



AdS/j^ soliton 



For completeness of the discussion on AdS solitons, we also give the results of the simplest 
d = 3 AdS soliton. The UV expansion of the solution r{z) is simply 



r{z) =R + a3{R)z^ + 0{z^ 



(25) 



where 03(7?) encodes the IR information and should be solved from the full equation of 
motion. From the above expansion we obtain 



dA 
dR 



-3as{R) + 0{e) 



(26) 



There is no UV-dependent divergent term in (26), but substituting (25) into the action ([T]) 



yields the divergent term 1/e. It is independent of R, hence does not appear in (26). 

For d = 3 AdS soliton, the cylinder solution is trivally r{z) = R, as could be seen from 
the equation of motion ([s]). As we will see, it is the dominated topology for large R. From 
the action (7) we obtain the on-shell action S = - — -. For this case we have 03 (i?) = 
and S'finite = ~^- 

We now concentrate on the disk solution which dominates over cylinder one at small R. 

The a^i^R) and Sgnitc for AdS4 soliton with f{z) = 1 — z^ is calculated numerically and 

plotted in Fig.[9} respectively. There is a phase transition at i? = 0.3686. For R < 0.3686 

the disk topology dominates, while for R > 0.3686 the cylinder topology is dominant. 

Based on the above result, we now apply the differential scheme (111) with d = 3 to 
^5(3) 
obtain — "^""^ for AdS4 soliton with f{z) = 1 — z^, and employ gsiRdajSa^ite to get the 

(3) I I 

corresponding 5'{jv_jj^j. Both of the results are shown in Fig. 10 Compared to the cases of 
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FIG. 9. Left: The as{R) for AdS4 soliton with f[z) = 1 — z^ . Right: Finite part of on-shell action. 
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FIG. 10. Left: The 



dS, 



(3) 
UV-i 

dR 



for AdS4 soHton with /(z) = 1 — z^ . Right: The corresponding 



5, 



(3) 
UV-ind- 



AdSs^e solitons, there is no violation of C-theorem even there is a first-order phase transition 
at -R = 0.3686. The renormaUzed entanglement entropy becomes constant after the quantum 
critical point. This is consistent with the expectation for the entanglement entropy of a 
(1 + l)-dimensional gapped system. 

B. Extracting the topological entanglement entropy 



According to the study of the strongly coupled condensed matter systems, the entangle- 
ment entropy contains both the short-range and the long-range ones |15] - H^ . The short-range 
entanglement is responsible for the area law nature of the entanglement entropy which mea- 
sures the number of the entangled pairs with one particle inside the chosen region and the 
other one outside. On the other hand, the long-range entanglement is a constant topological 
invariant, which is independent of both the UV and IR scales, and should be associated with 
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existence of the topological order. Especially, there are some exactly solvable model with 
topological order in (2+l)-dimensions, and their entanglement entropies have the structure 

S = aR--^ (27) 

where a and 7 are some constants. A nonzero 7 encodes the quantum dimensions of the any- 
onic excitations in the topological ordered phase, and is called the topological entanglement 
entropy. See [19] for the discussion of the topological entanglement entropy for the higher 
dimensional theory, which again should be a constant piece in the entanglement entropy. 

Since the topological entanglement entropy should be independent of the UV and IR 
scales, it should be encoded in the UV-independent piece. Note that the differential operator 
for AdS solitons ([s]) indicates 

_\wwRMRm-^)---iRm-id-^))s. dodd, 

•JUV-ind — \ \^0) 

[wwARm + ^)(RjR - 1) ■ ■ ■ (^ Jj - (^ - 3)) S, d even, 

hence the topological term would survive the differential operator in even dimensions, while 
in odd dimensions we could never observe such term in S'uv-md- Let's focus on AdS solitons 
in even dimension below. 

(d) ds'-'^'' 

To obtain •S'uy^j^j by integrating — ^^""^ over R, one will get an integration constant. 
However, this constant is not relevant for topological order since it can be fixed by the UV 
part of the UV-independent piece, namely, S\jy_^^^{R = 0). To look for the topological 
entanglement entropy encoding long-range entanglement, one instead should look for the 
IR behavior of the UV-independent piece. More precisely, one should extract the constant 
piece in the large R expansion of ^uy.jjjj. This piece will be independent of both UV and 
IR scales and should encode topological order. 

To avoid the numerical uncertainty, we here introduce an analytic way to extract the 
topological piece of entanglement entropy. The method is to consider the large R expansion 
of both action (JTl) and the associated equations of motion, and then order by order solve r^'s 
which are the coefficient functions in 

r{z) = r,{z)R + n(.) + ^ + 0{^) (29) 

where r^'s satisfy the UV boundary condition ro(0) = 1 and rj^o(O) = so that r(0) = R. It 
is easy to see that rQ{z) cannot be nontrivial from the leading order of equations of motion. 
Thus we set Tq^z) = 1. 
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To be specific, we consider d = A case. The action ([T]) in the large R expansion is 

A . r i. f ^^72^ + -^^ + ^SrL= + 0(i)] . (30) 

and the equation of motion in the large R expansion yields 



Solving (31) with the boundary condition ri(0) = yields ri{z) = by using the fact that 



/(O) = 1 and ri(0) is finite such that ri(0) = 0. This implies that the i?-independent term in 



(30) is zero irrespective of the topology of the holographic entangling hypersurface because 



we only use the UV geometry to yield ri{z) = 0. From (28), we conclude that the topological 



entanglement entropy is zero for both extremal and non-extremal AdSs solitons. 

III. CONSIDERATIONS FOR THE ADS4 BLACK HOLES 

We now consider another Lorentz non-invariant setting by turning on the temperature 
and chemical potential for the dual CFT. This is just to consider the AdS black hole with 
the following metric [38] (to be specific we consider the AdS4 planar black hole) 

ds' = ^ (-f{z)dt' + ^ + dr' + r'dA (32) 



with 

m = 1 - (1 + #^)(-)^ + #^(-)^ (33) 



272 '^z+' 272 ^Z+' 



„2 _ e'L\^s 



where /i is the chemical potential for the dual CFT and the parameter 7 = — ^ is the 
dimensionless ratio of the Netwonian and Maxwell couplings. Moreover, the temperature 
T of the black hole or the dual CFT is related to the position of horizon z+ and chemical 
potential fi by 





^ = 7^(3- 


272^- 




(34) 


The extremal black hole has T 


= by choosing 


^2 ,,2 

272 '^■ 






The thermal entropy density of the dual CFT is given 


by the Bekenstein-Hawking 


area 


law. 












2n Ah 

■Sthermal o t a 

fi;2 V2 


^2^2 

K Z: 


i 


(35) 
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where V2 is the field theory volume and Ah is the event horizon area. 



Now consider the holographic entanglement entropy in background (32). It is given by 
the area of the minimal surface determined by the action 



A= I v/det^iM= / dz^Jj+r^:= I dz C. (36) 



From (36) we obtain 



11=^ = ^ , n = Ur-C = , "" (37) 

dr ,2 JV^2 z^Vm + fr') 



f 
which appear in ^. 

For the AdS black holes, one may also expect that there are minimal surfaces of either disk 
or cylinder topologies. Since the horizon is just the coordinate singularity, one might expect 
that the minimal surfaces of "cylinder topology" would extend into the region inside the 
horizon and have turning points there, that is, they are in fact of disk topology with the tip 
shadowed behind the horizon. However, it was pointed out in [39J that the minimal surfaces 
could not extend all the way to the horizon, hence they cannot penetrate the horizon. We 
therefore just focus on the solutions with disk topology. 

Solving the equation of motion for the minimal surface in the UV expansion, we have 

riz) = R-^ + «3(i?)^' + 0{z') (38) 

where the higher order terms are not relevant for UV-independent piece of the entanglement 
entropy, and again a^{R) should be obtained by solving the full range of the equation of 
motion, and depend on the IR behavior of the minimal surface. Using ([9]), we have 

d A 1 

-^ = --ma,{R) + 0{e^). (39) 

To extract the RG fiow behavior of the renormalized entanglement entropy, we can ap- 
ply the differential subtraction scheme given in p2j. By using the commutator relation 
[Or, Rdi^ = dji, we obtain 

7 q(3) BH j , 

~W^ = RdR^ = RdR i-SRa^iR)) . (40) 

We first consider the case of non-extremal black hole. We solve 03 (i?) numerically and 



the result is shown in the left plot of Fig.jTT] The on-shell actions with divergent part 
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FIG. 11. Left: The a^^R) for non-extreme AdS4 black hole with f{z) 
part of the on-shell actions. 



l-2z^ + z^. Right: Finite 



dSuv-ind/dR 
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10 



FIG. 12. The 



dS, 



{3)BH 
UV-ind 



dR 



for non-extreme AdS4 black hole with f{z) 



l-2z-^ + z'^. 



5*^;^ ~ R/e subtracted are shown in the right plot of Fig. Kill Since there are only 
solutions with disk topology, there would be no phase transition along the RG flow. 
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We see that 



dS, 



{3)BH 
UV-ind 

dR 



The RG flow of the renornialized entanglement entropy is then followed from (|40|) and 
the numerical result is shown in Fig. 



is always positive, implying 



that more and more states are thermally excited as we go to higher temperature regime, or 
equivalently lower energy scale. For large R it is consistent with the linear running of the 
thermal entropy, i.e.. 



dS, 



(3) 
thermal 



dR 



ATlrC Sthermal- 



(41) 



Hence there is no phase transition but a smooth crossover interpolating between the renor- 
malized entanglement entropy for the ground state in the IR regime and the thermal entropy 
in the UV regime. This fact supports the conjecture proposed in [35] Q 

^ Their conjecture refers to the finite part of the entanglement entropy, the ^finite- Instead, we are consid- 
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FIG. 13. Left: The a3{R) for extreme AdS4 black hole with f{z) = 1 - 4z^ + Sz"^. Right: The 



corresponding 



dS, 



(3)BH 
UV-inc 

dR 



We also consider the case of extremal black hole, its 03 (i?) and the RG flow of S", 



UV-ind 



are solved numerically and the results are plotted in Fig. 13 Again we see the crossover 
from the IR regime to the UV one. 



IV. CONSIDERATIONS FOR THE ADS5 SOLITONS AND BLACK HOLES 
WITH GAUSS-BONNET CORRECTION 

In this section we will consider the effect of the Gauss-Bonnet term to the refinement 
of the holographic entanglement entropy for both AdSs soliton and black hole. The bulk 
theory we consider is given by the action 

/ = - ,. ^„ / d^x,/gl-—+n + 



C 



GB 



(42) 



IQtiGn J -^y^ y L2 ■ - ■ 2 
where Xgb is the coupling constant for the Gauss-Bonnet term with the Lagrangian 

Cgb = n^up^n'^''"' - ^n^.n^'' + n\ (43) 

In the dual theory, the higher curvature terms correspond to some perturbation in the sub- 
leading order of inverse 't Hooft coupling. 

The Gauss-Bonnet coupling Xgb should be in the interval [0, 1/4] for the metric to be 
well-defined outside the horizon of the AdS black hole or the IR end-point of the AdS 
soliton. Moreover, for the black hole in (4 + l)-dimensional AdS-Einstein-Gauss-Bonnet 



ering the UV-independent piece, the 5'uv-ind- In this sense, our results yield a refined version of their 
conjecture. One main difference is that the area-law piece of S'finito in the UV regime is UV-ambiguous and 
cannot be included in iSuv-ind- In contrast, the volume-law piece in the IR regime is UV-unambiguous. 
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gravity theory, it was shown in [301 EI] that the dual CFT will violate microcausality and 
render inconsistency when Xgb > joo' ^^ ^^^^ then explore this effect to the refinement 
of entanglement entropy by studying — ^^'"'^ for various values of < Xgb < 1/4- At the 
same time, we will check whether the Gauss-Bonnet term would affect the crossover from 
the UV regime to the volume law in the IR regime. On the other hand, the boundary dual 
theory of the AdS soliton is non-relativistic so that the micro causality constraint may not be 
relevanlj^ we will simply pick a specific value of < Xgb < 1/4 in the following discussion. 
We will now first consider the RG flow of the renormalized entanglement entropy for the 
Gauss-Bonnet corrected AdS soliton, and then for the corrected black hole. 



A. Renormalized entanglement entropy for the Gauss-Bonnet corrected soliton 

The AdS soliton solution in (4+ l)-dimensional AdS-Einstein-Gauss-Bonnet gravity the- 
ory is given by the metric f\ [^ 

*' = i%;57i;) + I|^^^ + /W-j. '■ = 0.1.2 (44) 



where 




/W = ^X^|l-Wl-4Aci.<|l-(-) M' (45) 

/o = hm f{z) = ; , (46) 



^Ads = -jf^ , and 9 r^ 9 + Le, Lg = ttzq • (47) 

v/o 

For Xgb — ^ 0, the metric (44) reduce to the AdS soliton part of ([s]). Note that L is different 



from Z/AdS) and in the numerical study of this section we will set L = 1 instead of //Ads = 1- 



From (45) and (46) it is easy to see that Xgb should be in the interval [0, 1/4] so that the 



metric (44) has the well-defined Euclidean section for < ^ < 2:0- 



Consider a disk on the boundary with radius R, the induced metric of the minimal surface 
is given by 

dsL = L= ( Jj (£- r(^)= + -^) d^' + £- ^ #= + M ,,= j , (48) 



^ Though, there is Lieb-Robinson bound [331 IM] for the non-relativistic systems, which plays similar role 

as the speed of light constraint on the signal propagation for the rclativistic one. 
^ In |16| , the UV divergence structure of the holographic entanglement entropy of this metric for the stripe 

region has been studied. They also studied the entropic phase transition by varying Xqb- 
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where r and are the radial and angular coordinates of the disk respectively. The minimal 
surface is determined by specifying r[z). 



The holographic entanglement entropy in the Gauss-Bonnet gravity is given by minimiz- 
ing the functional [HI [29] 



A 



dxWh (1 + XgbL'TZ) + 2 XgbL' / dxWh /C , 

7 A •^d"/A 



(49) 



where TZ is the intrinsic curvature of the induced metric h; hj, is the induced metric on d'jA 
and /C is the trace of its extrinsic curvature. The second term is the so-called Gibbons- 



Hawking term. From (48) we obtain 






^AdS- 



L^ (2fr — zfr — 2zfr + z^fr 



WLl,s + L'f 



where 



q{z) = X 



GB- 



L [Afr — zfr — 2zfr 



(50) 



(51) 



Integrating the term q{z) on 7^ gives rise to a surface term which cancels the Gibbons- 



Hawking term in (49). Therefore, the functional we need to minimize is 



A 



dz 



dz C . 



L^r 



L 



AdS 



^Ads + L'fr' + Xgb 



L^ (2fr — zfr — 2zfr + z'^fr 



(52) 



The equation of motion for (52) turns out to be 



'.Z 



r / r ( 3 - XgbzJ) r + 4 Xgb 



+ L'Ll,,fr' 



-1 + 2 XgbI - 2 Xgbz/ + Xgbz^I) 

■f (-2zf2 + r(3f + 2zr)) 



+ 2f[z{-l + 2 Xgb zf)r' + r (r(-6 - 3 Xgb zf + Xgb z^f) 

^ . r • / , •_, 2 



+ z{l-2XGBzf)r 
^^ {r{z)(3r — zf) + zf 



+ L^LX^s [zff(2 r(l + Xgb zf) - 3 Xgb z'fr] 



(-2 + 3 Xgb zf + Xgb z^f) 



+ 4 Agb/^ {r{z){3f — zf) + zf{—f + 3 zf)) + 2/ izf (f{ 

■^XGBZ^f)r) +r(f{-3-6XGBzf + XGBZ^f) + z{l + XGBzf)f^^'\ . (53) 
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FIG. 14. Disk (blue), Cusp (Green) and cylinder (red) solutions of the minimal surface for AdS 
soliton with higher derivative correction with \gb = 0.05, zq = 1. Since there are infinite number 
of cusp solutions with the same i?, we plot a typical one. 

1. Solutions of the minimal surfaces 



First, we consider the solution of (53) in the small R regime. In this regime, the solution 



has a disk topology as the blue line in Fig. 14 Near r = it can be expanded as 



where z^ is defined as 



z{r) = Zrn + Z2r^ + 0{r^), 



Zm = z{r = 0). 



(54) 



(55) 



By plugging (54) in (53), we find that the coefficient Z2 satisfies the following quadratic 
equation, 

azl + bz2 + c = 0, (56) 

where 



a = -\1zlzl\lj,{-\ - 7 + 2\gb)\^zI\gb + 4(1 - i + 2AGij(e - 2))], (57) 

h = -2ztz^{l + 7)AGij[24AGij(-8 + S\gb + H) 

+ 4(-5 + 2A\gb - IQ\Ib + Se - 14AgbO], (58) 

c = S^A^.^ + 4z^Agb(13 - 20Agb - 70 + 3^ (l + ^\Ib " ^ + Agb(-5 + 3^)) , (59) 

and for simplicity, we introduce 7 and ^ as 



7 



^l-AX 



GB, 



(60) 
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FIG. 15. Relation between iP' — Aac and Zm with Xgb = 0.05 and zq = 1 
The discriminant D{zm) = b"^ — 4:ac of (56) as a function of Zm is shown in Fig. 15 There 



is Zm = Zd which satisfies D{zd) = 0. Therefore, the solutions of disk topology exist only for 



0<zm<zd i.e., 0<R<Rd^ R{zd) 



In the following we will take Xqe = 0.05 and zq = 1, and in this case Ra = 0.528. 



(62) 



Next, we consider the solution of (53) in the large R regime. In this regime, the solution 



has a cylinder topology as the red line in Fig. 14 Expand the solution near z = zq as 



where 



and 



For 



r{z) =ro + ri{z - Zo) + O {{z - Zq)') 



'^o = rizo), 



ri 



-ro7' + 7 V^o7' - 3^o'(l + 7)Agb(1 + 8Xgb) 



UzqX 



GB 



v/3^o(1 + 7)Agb(1 + 8Agb) 

7 



(63) 

(64) 
(65) 
(66) 



the expression inside the square root in (65) becomes zero. Thus, the solutions of cylinder 
topology exist only for 

r-cyi < To i.e., Rcyi = i?(rcyi) < R. (67) 

In the case of Xqb = 0.05 and Zq = 1, i?cyi = 0.963. 
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For Rd < R < Rcyi, solutions of (53) have a cusp shape as shown in Fig. 14 For this 
solution, z'{r = 0) 7^ 0. However, the cusp solutions for a fixed R are not unique because 
we can adjust both z{r = 0) and z'{r = 0) to have the same R at UV. The absence of 
the smooth solution and the non-uniqueness of the cusp solutions suggests that there is no 
saddle point for prescription of [lOj in evaluating the holographic entanglement entropy in 
this regime of R. This may suggest the need of some quantum version of prescription of [10] 
to deal with such a case. Since we do not have such a prescription yet, in the following we 
will just skip discussion of the RG behavior for this regime. 

2. Renormalized entanglement entropy and its RG flow 
Recall ([9]) for the RG flow of the on-shell action. 



dA dz^ dr{e) 

dR = -^^'-^lR-''^'^^R^ (^^) 



where 
or 



L' {z^\gbL\J + fLl.si-'^zXGBLlds + L'r{l + zXcbDv) + L^ prr{-2\GBL\ds + LH')) 

(69) 



and 



n = Iir-C 



L^ {LH\GBf{-2f + zf)r^ + r(Li,5(l + 2Agb/ - zXaBf) + L^f{l + 4Agb/ - 2z\GBf)r 



rp^ 



z%L\,, + L^fT^f/^ 



(70) 



After simplification, the first term in (68) becomes 



'^(2:m)-7^ = , nn — N ^' ^°^ ^^^^ topology, (71) 

'H{zm)—r^ = 0, for cylinder topology because —^ = —— = 0. (72) 

Note that it is not zero for the disk topology, unlike the case with Xgb = 0. 
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The UV behavior of the solution r(z) is 



r{z) = R + 02-2 + a4^z + b^z log — H , (73) 

K 



where 



, ^ -1 - 7 + (7 + 57)AcB - 2(7 + 37)A^^ + (7 + 7)A^^ 

32i?3(-l-7 + (3 + 7)AGB) ■ ^ ^ 

Again the coefficient 04 cannot be determined from the UV expansion, and should be solved 

from the full equation of motion. 



Plugging (73) into (68), we obtain 



dA dz 

= --H^zn,)-^ + KRai{R) + UV-dependent terms + C(e), (76) 

dR dR 

where 

_ 3 256V27(-1 - 7 + (7 + 57)Agb - 2(7 + 37)A^^ + (7 + 7)A|^) 

(l+^)13/2(_l_^+(3 + ^)^^^) • I ) 

and 0{e) terms vanish at e — t- limit and are not relevant. The UV-dependent terms are 



a4 
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FIG. 16. The a4{R) of disk (blue)topology solutions for < i? < Rd = 0.528 and cylinder (red) 
topology solutions for Rcyi = 0.968 < R with Xgb = 0.05, zq = 1. For R^ < R < Rcyi, the solutions 
have cusp shape but are not unique. It suggests the absence of saddle point. We thus leave it open 
in the plot. 

29 



where 



and 



Cl 



_^3^(-l - 7 + (3 + i)Xgb + 4(1 + ^)XIb) 



(1 +7)5/2(1 + ^_4Agb) 



(79) 



C2 



-L' 



8^27^ 



(1 + 7)^2(1 + ^ _ 4Agb)=^(-1 - 7 + (3 + i)\gb) 
X [-1 - 7 + (11 + 97)Agb - 4(11 + 7i)\Ib + 7(11 + 57)A^b 

-5(11 + 37)AS^ + (11 + 7)A^5]. 



(80) 
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FIG. 17. 'H{zm)^ with Xgb = 0.05, zq = 1 for disk topology solutions 



We solve 04(^2) and T-Li^z^ 



' dR 



(for disk topology only) numerically, and the results are 



given in Fig. 16 and Fig. [TTj respectively. On the other hand, for the cylinder topology, 

TJ ( ■y \dZrn 



0. Using the above numerical data, we can then apply the same differential 



subtraction scheme given in (18) to extract from (76) the RG flow of the renormalized 



entanglement entropy in this case, and the numerical result is shown in Fig. 18 



Compared with Fig. |4] for the Xgb = case, we find that the Gauss-Bonnet interaction 
changes the UV behavior of — '^^'"'^ ? bringing serious violation of the C-theorem in the UV 



regime. This should be due to the contribution from T-i{zm)^^ in (76), which is dominant 



in the UV regime. Moreover, the Gauss-Bonnet term introduces the cusp solutions to make 
the transition between UV and IR regimes ill-defined, at least numerically. Despite that, 
the feature that the renormalized entanglement entropy tends to constant in the IR regime 
is still retained. 
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FIG. 18. Left: The 



dSu^ 



dR 



for the disk (blue) topology for < R < Rd = 0.528 and the cylinder 



(red) topology for i?cyi = 0.968 < R with \gb = 0.05, zq = 1. For Rd < R < Rcyi, solutions 
have cusp shape. In this region, there is no way to find the unique — ^Jr"^'^ with fixed R. Right: 
Zoom-in of the region around the phase transition. 

3. Extracting the topological entanglement entropy 

In order to extract the topological entanglement entropy, we perform the large R expan- 



sion for the action (|52|) and equation of motion (|53|) as following: 
A = 



^^ J L\Ll,,il + 2Xanf - XoBzf) + L^frj) ^ 



^AdS ^^ V -^AdS + -^^ / ^1 



+ 



L' 



^AdS 



.3 ( 7-2 



(^ids + L'fr!)^/' L 



ri(Li,s + LVrf)(Li,s(l + 2Agb/ - XcBzf) + L'fri) 



+ h [LU Xgb z'f + L'f\L'rir2 - 2Li,s ^gb {zt{ + f^)) 

^L\^^j{-2L\^^\GBZ^L\T^^\GBZf{zT\^T^))))\ + (^(^) ^ (81) 

= {LV'r?(-6/ + z/) + L^LLs/r2(^12AGB/'ri + z/ri(3-AGB^/) 

+8 \gb z fh + 2 /(- 6 ri - 3 Xgb zfh+ zh-2 Xgb z^fh + Xgb z^ fh) 
+ 2 L^Lids ( zfh + Xgb z^f^h + 2 Xgb fi^h - zh) + f {-^h - Q Xgb z fh 

+ zr, + Xgb z'fh + Xgb z'fh )) } i? + 0{R^) . (82) 
For the cylinder topology which dominates at large R, we have ri finite when 2 — )■ , for 



which equation (82) gives 

n(0) {[fohior + (1 - XgbIo)]" - AL/o} = . 
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(83) 



Since the term in the curly braces of (83) are positive definite, we have ri(0) = 0. Note that 



ri(0) = and the fact that ri = is a solution of equation (82). We then conclude that the 



unique solution to (82) is ri = . Then, it is straightforward to see that the i?-independent 



terms in (81) vanish. This yields zero topological entanglement entropy. This is consistent 

that the topological order will not show up in the leading order 



with the expectation in 

of 1/A^ expansion, which captures up only classical phenomena and not the quantum ones 

such as the topological order. 



B. Renormalized entanglement entropy for the Gauss-Bonnet corrected black hole 



Now we turn to case of AdSs black hole with Gauss-Bonnet correction. The bulk theory 
is the same as for the Gauss-Bonnet corrected AdSs soliton, and the black hole metric is the 



doubled Wick rotation of (44), which takes the form as 
L2 



ds' 



1 



, . fiz)dt' + -j—^dz' + /o(rfr^ + r\de' + sin^ 
z^ \ f{z) 



(84) 



where, f{z) and /o are the same as (46). 

By considering a disk on the boundary with radius R, the induced metric of the minimal 
surface becomes 



ds 



ind 



L' 



+ Ur{zf ) dz' + M^(rf^2 ^ gi^2 ...2^ 



^^2 \f{z) J Z 

where r, 6 and (f) are radial, polar and azimuth coordinates respectively. 



(85) 



The functional for the entanglement entropy is the same as (|49|). Using (|85|) we obtain 



V/i (1 + XgbL TZ) = V^^^ \- ^GB 



r^ , i^z"^ + fof{r'^ -2zrr + 2z^P)) 



z'^ 



z- 



WVl + hfr^ 



where 



q{z) = X 



GB- 



^fo^ r{r - zr) 
zW^ + fofr' 



+ g, (86) 



(87) 



Integrating the term q{z) in the bulk yields a surface term cancelling the Gibbons- Hawking 



term, then (49|) becomes 
A = 



, , 2/orVl + /o/r2 4{z^ + f,f{r^ - 2zrr + 2z^r^)) 
dz I — ^= h agb 



'Vf 



^y/v^l + fofr' 



dz C, 



(88) 
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from which we derive the equation of motion for r{z) as follows 
1 







foi-Qz'XcBizfr - 2fofr' - 2/(r - zr)) + Azr{l + zXgbJ 



+ ^^GBflfr^ - 2/(Agb + /o(-l - zXGBfY^) + kPr{2\GBr + hr^ - QzXgbt)) 
+ for\-zfr + /(6(1 + zXGBf)r - zfofr^ - 2zr) + 2fofr\-6XGBr + S/or^ - AzXgbt) 

- 2f{6XGBr - 6for' - 2zXGBr + zfj^r))). (89) 



The UV behavior of the solution riz) is obtained as 



r{z) = R + 



Agb(-1-7 + 4Agb) 2 , /u\ i 



2R{-1 + -f + 4Xgb] 



-z' + a^{R)z^ + ■ ■ ■ 



(90) 



where 7 is defined in (60) and a^iK) should be determined by solving the full equation of 



motion (89) 



Unlike the complication for the Gauss-Bonnet AdSs soliton case, there are well-defined 
solutions of disk topology for all R. As for the AdS4 black hole, we just need to consider 



the disk topology. The RG flow of the on-shell action is given by (68), and we need to see 



if the first term in mm has no zero contribution or not. From (88) we can obtain 



n:= 



5£ _ 2/ox/7(-4^AGBr + 2z''XGBf{Z + 2/0/7-^) -f hr''r\\ + ]{-2Xgb + /or'))) 



br 



z\\ + /o/r2)3/2 



(91) 



and 



'H:=nr-i: 



For disk topology. 



2{2z''Xgb - ^zXGBflPrr^ + kr\^ + ^^gbUPt^ + /(2Agb + /or'))) 



dr 

dz 



By plugging ( 93 ) into ( 92 ) , we get 

T-L{Zri 



00, 



\GB 



'v^(l + /o/r')3/2 



r{zm) = 0. 



0. 



forizr, 



(92) 
(93) 

(94) 



Then the first term of (68) becomes zero 



From (68), (90) and (91), we get 



dA 



dR RV2-fXl 



GB 



7 



A 



GB 



[-A^s (1 + 7 - 2(4 + 37)Agb + IQXls) 



-4R^-1 + 7 + 2Agb)(-1 + 4XGB)a^{R)] 
+ UV-dependent terms + 0{e) 



(95) 
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where the UV-dependent divergent terms becomes 

R 



27(1 + 7 + 4(-1+7)Agb) 1 



(96) 



Xgb (1 + 7-4Agb) e2' 

In fact in the action there is an additional R-independent logarithmic UV divergent term, 



which does not appear in (95). To see this, we substitute the UV expansion of r{z) (90) into 



the action (88) and find it as 



^^^ (1 + 7 - 12\gb) log - 
2 - 27 ^ ' ^ ' ^ R 



(97) 



which should be subtracted along with the quadratic divergence when evaluating 5'finitc- 
We then numerically solve 04 (i?) and the finite part of the on-shell action ^finite for differ- 



ent values of < Xgb < 1/4, and the results are plotted in Fig. 19 and Fig. 20 , respectively. 



dSuV-ind 



To extract the RG flow — '^^'"'^ from (|95|), we again apply the differential subtraction scheme 
based on [22] • Explicitly, it is 



dS- 



UV-ind 

~dR 



1 rIA 

-(i?9H + l)(i?9^-l)^. 



(98) 



The numerical results are shown in Fig. 21 

The Gauss-Bonnet interaction corresponds to some operator at the sub-leading order in 
the inverse 't Hooft coupling expansion. It denotes the finite coupling correction to the 
infinite 't Hooft coupling limit in the dual field theory. First, we notice that for Xgb = 0, 
i.e., corresponds to the AdSs black hole of Einstein gravity, the UV behavior of — ^Jr''"^ ^^ 
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FIG. 19. The a4(i?)'s for Gauss-Bonnet corrected AdSs black holes with zq = 1 and Xgb = 
(purple), 0.05 (blue), 0.09 (cyan), 0.2 (green), respectively. 
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FIG. 20. The S'snitc's for Gauss-Bonnet corrected AdSs black holes with zq = 1 and Xgb = 
(purple), 0.05 (blue), 0.09 (cyan), 0.2 (green), respectively. 
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FIG. 21. Left: The — ^^'"'^ 's for Gauss-Bonnet corrected AdSs black holes with zq = I and Xgb = 
(purple), 0.05 (blue), 0.09 (cyan), 0.2 (green), respectively. Right: Zoom-in of the region with 
the sign-change. 



different from the one for the AdS4 black liole case. In tlie latter case, the — ^^'"'^ is always 
positive, but liere it is negative at UV and tlien turn to become positive to capture tlie 
volume law nature at IR. We can see that the effect of tlie the Gauss-Bonnet interaction is 
to cliange tlie UV behavior of — ^^""^ so that it becomes all positive even at UV for large 
enough Xgb- 



Despite the discrepancies in the UV behaviors for different Ag^'s, the qualitative behav- 
iors of the results are the same: the on-shell actions catch the volume law of the thermal 
entropy in a smooth way as R becomes large, just like the AdS4 black hole case. This again 
supports the postulate proposed in [35j . We thus conclude that the crossover is not effected 
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by turning on the Gauss-Bonnet interaction. 

In [Sni EI] it is pointed out that the holographic dual field theory with Xgb > 
violate microscopic causality, however, although the small R behaviors of S'finite and RG 



^ will 



fiow become quite different for sufficiently large Xgb, e.g., Xgb = 0.2 in Fig. 20 and Fig. 21 



nothing exotic happens in this regime. This agrees with the same consideration for the 
Gauss-Bonnet corrected AdSs soliton in pSj. However, there are some concern about the 
relation between the quantum entanglement and the causality formulated from the con- 
sideration of the quantum information sciences [HI |12], it may deserve further study to 
understand this issue in the context of holographic entanglement entropy. 
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FIG. 22. Left: The a2{R) for Xgb = 1/4 and zq = 1. Right: The corresponding RG flow of the 
refinement. 



Finally, we would like to give the numerical results for the Xgb = 1/4 case, for which the 
viscosity to entropy ratio vanishes for the holographic dual field theory. The peculiar feature 
of the geometry is the harmonic function f{z) becomes 1 — (— )^, which is quite different 

' — Y. The UV behavior of the solution is 



from the one for asymptotic AdSs, namely, 1 
different from the Xgb < 1/4 cases and is given by 



r{z) = R+a2Z^^ 



2 a2Z^{3{z^ + IGa^z^) + 64i?(a2^^ + 12a^4) + 4i?2(l + 96ai4 + 768a|4)) 



4^2(3^2 + 4Sa2Rzi + 4i?2(l + 480^^2)) 



(99) 
where a2{R) instead of a^{R) should be determined by solving the full equation of motion. 



and the result is shown in the left plot of Fig. 22 From (99) and (68) we can obtain 



dR 



2lV2{{4 



)a2 + ASRal + 128R 



2^3 



rA 



2^))]+0{e^). (100) 



The result of the RG flow is shown in the right plot of Fig. 22 , from which we see again the 
volume law for large R as well as the crossover. 
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V. CONCLUSIONS AND DISCUSSIONS: IR FIXED-POINT STATE FROM ADS/MERA 

In this paper, we have considered the refinement of the holographic entanglement entropy 
and its RG flow behavior for the the non-relativistic systems of AdS solitons and black holes. 
The holographic entanglement entropy has different scaling behavior from the relativistic 
cases, so does the UV-independent piece, i.e., the renormalized entanglement entropy. As 
for the cases of AdS black holes, our results yield the refined version of the conjecture 
given in [35j that the transition of the UV cutoff-independent piece of the entanglement 
entropy between the IR and UV regimes is a smooth crossover even with the correction of 
the Gauss-Bonnet interaction. On the other hand, for the AdS soliton cases, we find that 
the renormalized entanglement entropy for ci = 4, 5 is not monotonically decreasing along 
the RG fiow, suggesting the violation of the C-theorem around the phase transition. Nor are 
they always positive definite. This indicates that the renormalized entanglement entropy 
might not be a good measure of the number of degree of freedom for d > 3. 

Turning on the Gauss-Bonnet interaction will only affect the UV but not the IR behaviors 
though the region around the confinement phase transition becomes ill-defined. Similarly, 
the irrelevance of the Gauss-Bonnet interaction to the topological entanglement entropy is 
also checked so that there is no non-trivial topological order for the AdSs soliton and its 
Gauss-Bonnet cousins. 

Before ending this paper, we would like to devote the rest of the discussions on how 
to understand the entangling nature of the IR fixed-point states of the holographic dual 
theory based on the conjecture of AdS/MERA proposed in [501 El]- We will argue that 
non-extremal AdS soliton has the product state as its IR fixed-point state, and the extremal 
AdS soliton instead has the nontrivial entangled state as the IR fixed-point state. The 
different nature of the IR fixed-point states depends on the topology of the large R entangling 
hypersurfaces. If our arguments here hold, this may be seen as another triumph of AdS/CFT 
in using the simple geometric picture to characterize the entangled mean field states. Further 
development along this line may reveal the holographic and geometric classification of the 
topologically ordered phases in the strongly interacting condensed matter systems. 

Though the wave function of a many-body system could look quite complicated, it could 
be simplified a lot through some appropriate local unitary operations, especially when these 
operations are adopted to remove short-range entanglement among neighboring particles. An 
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example of such unitary operations is the CZ (controlled-Z) operation, which transforms a 
Bell state into product state as 

CZm\+) + |1)|-)) = (|0) + |1))|+) = v^H-)|+), (101) 

where |±) = 4=(|0) ± |1)). Moreover, if we are only interested in the low energy behaviors of 
the system, we could further coarse-grain the wave function by merging the the neighboring 
sites after removing the short-range entanglement. After repeating the above two steps, we 
will obtain a far more simple wave function at the IR fixed-point, or the so-called mean field 
state. This is the so-called quantum state RG transformation [SH 1^3] (see also [Sni EZ] for 
practical numerical study) as shown in Fig. [II and can be adopted to classify the phases 
of the many-body systems. That is, all the wave functions flowing to the same fixed-point 
state under quantum state RG transformation describe the same phase. According to this 
scheme of classification, for gapped systems one may expect two kinds of the IR fixed-point 
states. One is the product state which encodes no quantum entanglement. The other kind 
is the nontrivial topological ordered states, which encode either long-range entanglement or 
some short-range entanglement protected by symmetries [2S] . In this way, one can tell which 
phase the system belongs to by looking into the IR fixed-point wave function, instead of the 
UV ones. In other words, the gapped systems are classified by the patterns of the quantum 
entanglement of the IR fixed-point states. Especially, for 1-dimensional spin chain, it was 
shown that all the ground states will flow to trivial product state under generic quantum 
state RG transformation unless some symmetries are preserved during the RG flow [^31 - 
HH]. However, the classiflcation of higher dimensional systems are still under development. 
The above scheme of looking into the IR fixed-point state is in contrast to what has been 
adopted in this paper and summarized in ([I]) by looking into the UV scaling behaviors of 
the entanglement entropy for the relativistic CFTs. 

The local unitary operation and the coarse-graining in the quantum state RG trans- 
formation can be implemented as the quantum gates of the quantum circuit with some 
pre-prepared inputs. Therefore, the whole procedure can be viewed as some time evolving 
procedure and then be implemented to solve some many-body systems. This idea then 
results in algorithm of multi-scale entanglement renormalization ansatz (MERA) [52] , and 
see [53] for more detailed introduction. In MERA, the local unitary operations in remov- 
ing the short-range entanglement are called disent anglers, and the merging operations for 
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FIG. 23. MERA network and its dual AdS geometry. Here the disentanglers are denoted by solid 
squares, and the isometries by the solid triangles. The links at different levels encode short-range 
entanglement at different scales. (I) The MERA for GET and its dual is the AdS space. Note that 
the depth of the MERA can be extended indefinitely as its dual AdS geometry. (II) MERA for 
gapped system and it dual is the AdS soliton. The MERA and its dual geometry end at some IR 
scale. Eor simplicity, we just plot the one-dimensional MERA, however, it is straightforward to 
plot for higher dimensional cases. 

coarse-graining are called isometries. Then, the whole procedure of quantum state RG 
transformation can be piled up as a network of disentanglers or isometries. The depth of 
the MERA network can be thought as the time evolution or RG flow, and the links in the 
network denote the short-range entanglement among the neighboring sites. A typical MERA 



network for both CFT and gapped system are depicted in Fig. |23j Note that the depth for 
the CFT is indefinite due to the scaling invariance and could be infinite for an infinite UV 
system. On the other hand, the depth for the gapped system is finite as the RG procedure 
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must end when reaching the IR mass gap. 

In practical, the MERA can be used to solve the ground state of the system by treating 
the disentanglers and isometrics as the variational ansatz, which can then be determined by 
minimizing the expectation values of the Hamiltonian. For examples, see [53] for this kind 
of applications. 

On the other hand, the MERA network yields a geometric picture of the quantum state 
RG, and indeed the geometry can be characterized by the aspect ratio of depth to width, i.e., 
z ~ log \x\. This aspect ratio encodes the block decimation of coarse-graining and is roughly 
coincident with the AdS geometry as first observed in jSH] and made more precise later in 
[51] . For the gapped system, the finite depth is consistent with the geometry of AdS soliton 
with zq ~ log^ where ^ is the correlation length. Moreover, by utilizing the unitarity feature 
of disentanglers and isometrics in the MERA network one finds that a site is only affected 
by the sites within its causal cone. The correlation between two distant sites are encoded 
by the intersection of the causal cones, which is pretty much the same as the geodesic in the 
AdS bulk connecting two boundary points. This then reminds the prescription of evaluating 
the boundary correlation functions in the AdS/CFT correspondence [13]. By the aspect 
ratio of depth to width, the length of the intersecting causal cone then yields the expected 
power law for GET correlation function and the exponential decay behavior for the gapped 
one. 

Similarly, the geometric picture of the holographic entanglement entropy is encoded in 
the minimal surface covering the boundary sites inside the chosen region as depicted in 
Eig. 24 for both GET and gapped systems (see also Eig. fl]). The entanglement entropy 



is proportional to the number of the links intersecting with the minimal surface because 
the links carry the short-range entanglement between the sites just inside and outside the 
chosen region. This then results in the expected area law for both GET[_]and gapped systems. 
The most interesting point is that the link at different depth level of the MERA network 
actually encodes the short-range entanglement at the corresponding scale. To be more 
specific, the links at level (the UV boundary) encode the short-range entanglement between 
nearest neighboring sites, but the links at level 1 encode the short-range entanglement 
between the next-nearest neighboring sites with the distance measured by the UV scale. 
Therefore, MERA network geometrically and systematically displays how the short-range 
^ It can also recover the logarithmic behavior for the 1 + 1 CFT. 
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FIG. 24. Minimal surfaces for the entanglement entropy in the MERA and in its dual geometry. 
The entanglement entropy is obtained by counting the links which intersect the minimal surface. 
This implies that the entanglement entropy is contributed by the short-range entanglement at all 
length scales smaller than the linear size R of the chosen region. (Ill) For the CFT case, the 
minimal surface is always in the disk topology. (IV) For the gapped systems dual to non-extremal 
AdS soliton, the topology of the minimal surface changes from the disk at smaller R to the cylinder 
at large R. Compare the minimal surfaces for MERA and AdS soliton, we conclude that the IR 
fixed-point state (the yellow part excluded from the minimal surface) is a product state since the 
links at the top level of MERA have no counterparts at the geometry side. On the other hand, 
for the extremal AdS soliton, the minimal surface is always in disk topology, this is because the 
extremal AdS soliton has gapless KK modes which may retain the entangled pairs at the IR regime. 
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entanglements of different lengtli scales are contributed to the total entanglement entropy 
of a chosen region at UV level. 

Especially, for the gapped system there exists a top layer in the MERA network, which 
represents the IR fixed-point and also encodes the short-range entanglement of the IR fixed- 
point state. Therefore, if the fixed-point state is not a product state, its short-range entan- 
glement will contribute to the total entanglement entropy. Otherwise, there is nothing to 
contribute for a product state. This then corresponds to the following geometric picture. 
Due to the existence of the IR top layer, the minimal surface covering the chosen region 
will have a fiat bend-over near the top layer. If the fixed-point state is the product state, 
then the fiat bend-over region of the minimal surface collect no entanglement from the fixed- 
point state. In this case, fiat bend-over region can be effectively removed, and the resultant 
minimal surface can be effectively viewed as ending on the top-layer. This is indeed the IR 
dominating cylinder topology found in the non-extremal AdS soliton case. From our above 
argument, it implies that the IR fixed-point state is the product state. This result is con- 
sistent with the vanishing topological entanglement entrop}|j and the negative value of the 
finite part of the holographic entanglement entropy, which could compensate the positive 
UV contribution to make zero total entanglement entropy near IR fixed-point. 

On the other hand, for the extremal AdS soliton case we see that only disk topology 
exists so that the flat bend-over region does contribute to the holographic entanglement 
entropy. From the above argument, this could imply that the IR fixed-point state may 
not be the product state but a nontrivial entangled state. Geometrically, the difference 
between extremal and non-extremal AdS soliton is that the spectator U(l) cycle for the 
former becomes non-compact at the IR fixed-point. That is, the IR fixed-point is a gapless 
state as the KK modes become massless at the extremal limit. These gapless excitations 
may retain some entangled pairs at the IR scale indicated by the flat bend-over region. 

The above speculation of the entangled properties of the IR fixed-point state from 
AdS/MERA can be further exemplified by our study of the AdS black hole. In this case, 
the AdS geometry provides more useful information than MERA, whose finite temperature 
version is barely studied. Based on AdS/MERA, the finite temperature MERA network of 
the GET is no longer extended indefinitely but will be terminated by the IR scale fixed by 
® We restrict our discussions here for the AdS^ sohton case, which is dual to the 2+1 gapped system. On 
the other hand, the nature of the topological entanglement in higher dimensional system is not clear. 
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the temperature. This will be a helpful guideline when implementing the finite temperature 
MERA for CFT. Moreover, from our numerical study we see that the the dominant topol- 
ogy at large R is the disk one whose refined holographic entanglement entropy captures the 
volume law of the thermal entropy. According to the same consideration as for the AdS 
soliton case, this implies that the IR fixed-point state has nontrivial entanglement at IR 
scale. Indeed, the IR fixed point state should be a thermally mixed state and is different 
from the product state dual to the cylinder topology. Though we may need the pattern of 
thermal MERA to understand the how the multi-scale entanglements distribute at nonzero 
temperature. 
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